A hidden confining world on the 750 GeV diphoton excess 
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We explain the recent diphoton excesses around 750 GeV by both ATLAS and CMS as a singlet 
scalar ® which couples to SM gluon and neutral gauge bosons only through higher dimensional 
operators. A natural explanation is that ® is a pseudo-Nambu-Goldstone boson (pNGB) which 
receives parity violation through anomaly if there exists a hidden strong dynamics. The singlet and 
other light pNGBs will decay into two SM gauge bosons and even serves as the meta-stable coloured 
states which can be probed in the future. By accurately measuring their relative decay and the 
total production rate in the future, we will learn the underlying strong dynamics parameter. The 
lightest baryon in this confining theory could serve as a viable dark matter candidate. 


I. INTRODUCTION 


Very recently, ATLAS and CMS collaborations 
have announced their first hosts of new results 
based on 3.2 fb~' and 2.6 fb‘ integrated lumi- 
nosity at LHC Run II ys = 13 TeV [1, 2]. Among 
various channels in searching for new physics, 
there is an intriguing existence of diphoton excess 
around 750 GeV, with a local significance of 3.6 o 
and 2.6 ø respectively in ATLAS and CMS. With 
more data accumulating, whether this is due to 
a statistical fluctuation or some manifestation of 
new physics would be revealed soon. Neverthe- 
less, as theorists, we should always be aware if 
this could be the first light that changes our cur- 
rent understanding of microscopic physics. 


The first appearing of this anomalous dipho- 
ton resonance at LHC Run II would unambigu- 
ously tell us some information. First, due to 
the Landau- Yang theorem [3], this resonance can 
only be spin zero or two instead of one. Second, 
the resonance decay into diphoton process can 
only be through the higher dimensional opera- 
tors [4]. Therefore, an unsuppressed total decay 
width would require an unconventional large pro- 
duction rate and one might need to try hard to 
hide its main decay channel into the SM back- 
grounds. Third, according to the 8 TeV LHC 
Run I results, CMS search [5] sets a 95% CL ob- 
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served upper limit of o(pp > ®)Br(® > yy) < 
1.5 fb, and ATLAS search [6] also imposes a sim- 
ilar constraints on RS gravitons. In order to ac- 
commodate both LHC Run I and Run II results, 
a larger enhancement on the diphoton signal from 
8 TeV to 13 TeV is needed and the gluon ini- 
tial state is preferred. Collecting all the above 
hints, we consider a singlet scalar ® with only 
SM higher dimensional couplings to gluon and 
neutral gauge bosons as perhaps the most opti- 
mal solution. 


While the process gg => ® — yy looks sim- 
ple, it does have a very rich and deep physics 
behind it. If the ® is pseudo-scalar, or even a 
pseudo-Nambu-Goldstone boson (pNGB), then 
the only existence of higher dimensional cou- 
plings to gluon and neutral gauge bosons is a 
natural consequence of ® parity violation due 
to anomaly [7]. The anomaly induced process 
at the IR, which affects the ® production and 
decay, is proportional to the number of colour 
Ny in the underlying confining strong dynam- 
ics. Therefore if this excess continues to exist in 
the future, by accurately measuring the dipho- 
ton resonance rate and the relative rate among 
different SM diboson decay channels, we could 
learn N, and the hypercharge of the confining 
vector fermions. This provides us another ex- 
ample of learning the ultraviolet physics at the 
infrared just like a rediscovery of colour Ne = 3 
in QCD through 7 — yy. If one of the confin- 
ing vector fermions w is a SM singlet, the baryon 
which is made of N, copies of Y could be a com- 
posite dark matter candidate if w is the light- 


est confining vector fermions. Therefore, we use 
this particular choice of charge assignment as the 
benchmark of our model in the phenomenology 
discussion. 

The layout of this paper is as follows. In 
Sec. II, we parametrize our setup in terms of 
effective theory and give a numerical fit to the 
13 TeV diphoton excess. In Sec. II, we build 
up the simplest hidden QCD model with the re- 
quired singlet couplings through anomaly based 
on SU(4) flavor symmetry and discuss the related 
phenomenology. The benchmark model with a 
composite dark matter (DM) candidate from the 
hidden baryon is highlighted. We make conclu- 
sions in Sec. IV. We provide an Appendix A to 
derive the effective Wess-Zumino-Witten (WZW) 
terms of the simplest confining model. 


I. THE EFFECTIVE THEORIES 


We first consider the general dimension-five 
couplings between a singlet (pseudo)scalar and 
the gauge fields, 
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respectively. Here, Gi, and F),, are SU(3). and 
U(1)em field strength tensors. F is expect to be 
the energy scale of the underlying new physics 
model. The decay widths of (pseudo)scalar 
®5(®p) to gluon and photon pairs could be ob- 
tained from Eq. (1), 
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There can be different origins for obtaining 
the operators in Eq. (1). The first example 
is analogous to the Higgs effective couplings 
to gluons and photons. One may consider a 
set of vector-like colour-charged and/or electric- 
charged fermions F; coupling with the singlet 
scalars such that 5°, Fy (Agp®+i7sAn®p) Fr . In- 
tegrating out the heavy fermions Fk, the Wilson 
coefficients in the effective Lagrangian Eq. (1b) 


are calculated to be, 
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with N.(k) and Q, being the SU(3)e colour de- 
grees of freedom and charges carried by Fẹ and 
we assume degenerate masses of Mp = F for 
simplicity. In the following, we will denote the 
couplings and charges of vector-like fermion col- 
lectively as Ap = {Ax , Ax} and Qr = {Q}. In 
addition, the effective coupling for the C P-odd 
operator can be induced by the anomaly of the 
chiral symmetry breaking. This is nothing but 
the effective WZW term [8, 9]. In the next sec- 
tion, a specific model construction will be given. 

Due to the couplings, we expect the following 
relations [75] 


(C,/F)}, (4a) 
(C,/F)?. (4b) 


The cross section of the production process of 
gg —> ®s/p should only depend on (C,/F)?, 
hence the following semi-analytic formula for 
the signal process could be expected, olgg > 
Dsyp > 1] ~ (Cg/F) (C2 f2)/(C2 fa + C2 fr). 
Here (fı, f2) are constants and one expect 
fi/fe x (as/a)2 ~ (10? — 10%) (notice that 
the charge Qk is absorbed into Cy. We find 
alag + s/p] = (6.2 fb) x O2( HEY)? 

We perform the numerical analysis of the 
diphoton excess by using the implementation 
of the dimension-five operators Eq. (1) in 
FeynRules [10] and generate events with Mad- 
Graph [11], interfaced with Pythia [12] and 
Delphes [13] for the parton shower, hadronization 
and the fast detector simulations. The analysis 
is conducted based on the CMS cuts in Ref. [2]. 
The diphoton events are reconstructed by select- 
ing photons such that pr(y) > 75 GeV, |n(7)| < 
2.5 and |n(y)| not within (1.4442,1.566). At 
least one photon should be in the barrel region, 

e., |n| < 1.4442. The diphoton invariant mass 
should be m,., > 230 GeV. For events with one 
photon in the end-cap region such that |n(y)| > 
1.566, we require Myy > 320 GeV. Furthermore, 
only the events with the diphoton invariant mass 
in the range of my, € (650 GeV ,800 GeV) are 
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selected. In order to account for the diphoton 
excess at LHC, the cross section of o(pp > ® > 
yy) at the LHC 13 TeV should be around 3 — 13 
fb, while the Run I constrains tell us that the 
cross section at the LHC 8 TeV should be less 
than 1.5 fb. Fig. 1 shows the contour plot on 
the (C,/C,,C,) plane for the o(pp > ® > 77) 
at the 8 TeV and 13 TeV LHC, from which we 
can infer that small C, < 3 requires a large ratio 
C/C; Z 8. Notice that in the case of degenerate 
coupling and charge for the vector-like fermions, 
Cy/Cqg = N-Q%,/4, therefore from the plot, we 
can see that our theory even remains to have 
small coupling Ar < 1 for Qr Z 3 when we fix 
Mr =1 TeV (C3 = 4AF Nf/3 and cP = 2ArAR Nf 
through Eq. (3)) 
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FIG. 1: The cross sections of o[pp > ® > y7] 
(unit: fb) on the (C,/C,,C,) plane at 8 TeV (in 
green dashed) and 14 TeV (in blue solid) with F = 
1 TeV. The pink region represents the region that 
can explain the LHC 13 TeV results while remains 
unconstrained from LHC Run-I bound. 


IHI. THE MODEL WITH NEW STRONG 
DYNAMICS 


Next, we turn to a specific model setup, by as- 
suming a new QCD-like strong sector with the 
gauge symmetry of SU(N,,). The new strong dy- 
namics possesses the properties of confinement 
and the asymptotic freedom. We denote the pion 
decay constant of new strong dynamics as fn, 
and the dynamical scale as A,. They are re- 
lated as An œ 4r fn/vV Nn by the large N scaling 


relation. Unlike the technicolour theories [14— 
16], the strong sector is not necessarily related 
to the electroweak symmetry breaking. As a 
result, we are free to consider the case where 
only gauge bosons of unbroken gauge symme- 
tries, namely, gluons and photons, can talk di- 
rectly to the new sector. We assume a set of 
vector-like fermions (7! ,w?) under the funda- 
mental representation of the new gauge group 
SU(N,). We are especially interested in the case 
that only the SU(3), and U(1)y fields have the 
anomaly in order to account for the diphoton ex- 
cess. This can be realized by assuming that the 
fundamental fermions (Y1 , Y?) in the new strong 
sector belong to the singlets of SM SU(2)L gauge 
group. In order to embed the SU(3). group and 
have the colour-singlet pNGBs, we consider the 
minimal case with Ny = 4 and gauge the sub- 
goup SU(3). x U(1)y of SU(4). The quantum 
numbers of (7! , Y?) are summarized in Table. I. 
More general discussions of the global symme- 
try breaking patterns in different representations 
under the new gauge symmetries can be found in 
Ref. [17]. The LHC phenomenology of the mod- 
els with vectorlike confinement were studied in 
Refs. [18, 19]. 
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TABLE I: The field content of the minimal model. 
VYi,2 EQ. 


All gauge anomalies are cancelled since the 
fermions (y!,7?) are vector-like under both 
SU(N,) and the SM gauge groups. The require- 
ment of the asymptotical freedom of the new 
gauge theory SU(N,,) bounds the number of fla- 
vor to be Ny < 11N,,/2, which is easily satisfy 
for Ny = 4 and N, > 3 in our model [76]. In ad- 
dition, the asymptotical freedom of QCD should 
be retained, which bounds the color degrees of 
freedom in the new strong sector as Np < 10. 

In the limit of vanishing SM gauge couplings, 
the strong sector possesses the global chiral sym- 
metry of SU(4)r x SU(4)R x U(1)a x U(1)s,, 
where U(1)p, denotes the baryon number sym- 
metry in the new strong sector. The axial sym- 
metry of U(1)a is broken by the instanton ef- 
fects and will be neglected henceforth. The 


confinement of the theory at a new scale A, 
will induce spontaneous chiral symmetry break- 
ing of SU(4)L x SU(4)R —> SU(4)y, which re- 
sults in 15 pNGBs II“. They can be decom- 
posed into representations of SM gauge goup 
SU(3). x SU(2) x U(1)y as follows 


15> (8 ; 1)o + (3 3 l)y, -y> + (3 5 l)y,-y, + (1 ; 1)o R 


(5) 


Below, we label the colour octet, triplet and sin- 
glet pNGBs as ®g, ®3, and ®1, respectively. The 
masses of pNGBs arise from the gauge-invariant 
mass terms of —mo(Y!y! + Y?y?) for vector- 
like fermions. Such mass terms provide the 
dominant mass source and the gauging of SM 
group only splits the colour octets ®s and sin- 
glets @,;. The meson masses are given by the 
Gell-Mann—Oakes—Renner relation [20] 
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The singlet remains the lightest as a result of the 
Witten theorem [21] and is a good candidate for 
the diphoton excess. Below, we focus on a bench- 
mark model with fy œ 2.5 TeV with Nn = 10 for 
a composite DM candidate in the spectrum. The 
corresponding vector-like fermion mass is found 
to be mg ~ 1 GeV to accommodate a 750 GeV 
singlet in the new strong dynamical sector. Note, 
the benchmark model given here is made by as- 
suming the unitarity bound to the DM mass is 
saturated. More generic choices of the pion de- 
cay constant of fy and the gauge symmetries 
can be expected. Correspondingly, the vector- 
like fermion mass range spanning from O(1) GeV 
to O(10) GeV is also natural. 
The pion-number violating interactions arise 
from the effective WZW term 
N, II ~ 
Lwow = ae ge ee , 
(7) 


where gg,c are the gauge couplings associated 
with the SM gauge field strength tensors F ae 
and the dual field strength tensor is defined as 
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formed over the SM indices. 


In principle, we can freely choose the hyper- 
charges of SU(3), triplet and singlet. The singlet 
couplings to the hypercharge fields lead to the 
possible diphoton signals, given that Yı # Yo. 
Note that topological interactions between the 
colour-triplet and the SM gauge bosons are for- 
bidden by the SU(3), symmetry. An interesting 
and special case is Y> = 0 where wy? is a SM 
singlet. In this case, the baryonic composites 
(Y) = et trade wt could be a dark 
matter candidate. Therefore, we will use Yo = 0 
as our benchmark model in the later discussions. 


According to the effective WZW terms in 
Eq. (7) for the singlet ®,, we expect its couplings 
to the (gg,77,yZ,ZZ) final states. Likewise, 
the colour octet ®g can couple to (gg, gy, gZ). 
The explicit expressions for the WZW effective 
terms are given in Eqs. (A4) and (A5), with the 
corresponding partial decay widths of @; and ®g 
given in Eqs. (A7) and (A8). In Fig.2, we dis- 
play the decay branching fractions of the singlet 
®, with the varying inputs of |Yi| by assum- 
ing that Y> = 0. The increasing |Y;| apparently 
leads to the enhancement of the decay modes of 
(17,92, ZZ). 
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FIG. 2: The decay branching ratios into various 
final states for the singlet 1. 


If the future LHC experiments discover both 
singlet ®,; and octet ®g, one may measure the 
signal rates to yy and dijets for ®1, and the sig- 
nal rates of y plus jets and dijets for ®g. Obvi- 
ously, these signal rates can be used to determine 
the hypercharges in the minimal model. From 


Eqs. (A7) and (A8), one has 
5 as P[®s > g7] 


Yr = ’ 
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Note from the meson mass spectrum in Eq. (6), 
the octets are typically O(1) TeV heavier than 
the singlets. 

The total production cross sections can be ob- 
tained by mapping the EFT parameter C, into 


the minimal model, which is Cy = 272 


Eqs. (7) and (A4). It is crucial to note the pro- 
duction of the colour singlet ®, is proportional to 
the number of colours in the new confining strong 
dynamics. Based on our numerical simulations 
by following the approaches described in the pre- 
vious section and the decay branching ratios dis- 
played in Fig. 2, we plot the signal predictions of 
o[pp > ®ı > yy] on the (fun ,|Y1|) plane in Fig. 3 
by fixing Y> = 0 for SU(10) hidden gauge sym- 
metry. The parameter region with large |Y1| and 
small fq inputs have been excluded by the AT- 
LAS searches for Zy [30]. Furthermore, the total 
decay width of ©, is found to be hardly larger 
than 0.1 GeV for parameter regions favored by 
the LHC 13 TeV and 8 TeV data sets. Therefore, 
the limits on cross section of the spin-0 resonance 
to diphoton in the narrow width hypothesis with 
Tx =0.1 GeV [5] is applicable to some extent. 
In addition to the diphoton signal predictions 
we explored above, the minimal model also pre- 
dicts several other experimental signatures. We 
note above that the colour singlet ®; decays also 
to gg. Correspondingly, one would envision the 
future dijet searches around the mass resonance 
of ~ 750 GeV. Expressed in terms of the effec- 
tive couplings defined in Eq. (1), the ratio be- 
tween the diphoton and gluon pair signals are 
determined by (C,/C,)? = 3,(¥?—-Y?)?. There- 
fore, the future observation and measurements of 
the dijet signals are not only useful for justifying 
the model, but also crucial for determining the 
hypercharge differences for the underlying model. 
The other experimental constraints related to 
the minimal model are the searches for the 
vector-like quarks w' and the colour triplet ©; 
in the spectrum. They can be pair produced and 
hadronize with quarks and gluons to form the 
R-hadrons. The R-hadron searches at the LHC 
thus place mass limits on the colour triplet, and 
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FIG. 3: The diphoton cross section of a[pp > ®; > 
yq] (unit: fb) on the (fm, |¥i|) plane for Yo = 0, with 
K = 1.6 [31] for the NLO QCD correction. The blue 
and green curves correspond to the 13 TeV and 8 
TeV predictions, respectively. The total decay width 
of ®; in unit of GeV is shown by purple contours. 
Cyan region are allowed region of Zy final states from 
CMS Run I, and magenta region is favored by 13 TeV 
diphoton excess and allowed by 8 TeV CMS diphoton 
limits. The colour number in the hidden strong sector 
is chosen as Nn = 10. A benchmark point marked 
by star correspond to fn ~ 2.5 TeV, |Yi| = 4/3, and 
the diphoton signal cross section of 5.4 fb. 


the colour triplet mass is bound above around 
845 GeV by LHC Run-I results [32, 33]. The 
LHC Run-II at 13 TeV would set the bound more 
stringent [34]. The analysis in Ref. [35] set an ex- 
clusion to charged stable particles to ~ 900 GeV 
for sufficiently long decay lengths of cr > 10 m. 

The colour triplet ®3 has hypercharge Yı — 
Y2, so that it can decay to SM lepton-quark pair 
through higher-dimensional operators for specific 
choices of Yı — Ya. Examples include 


i _ 
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(9) 
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For the benchmark in Fig. 3, one may look for 
the pair-produced lepton-quark signals of 772@ to 
search for this resonance. The current mass ex- 
clusion of ® is ~ 1 TeV from the ATLAS 8 
TeV searches [36]. For sufficiently small A, the 
metastable ©3 is expected to decay before the 
era of the big bang nucleosynthesis (BBN) [37— 
39] through Eq. (9). This is different from the 
case in Ref. [28] where small symmetry-breaking 
Yukawa couplings between the new fermions and 
the Higgs doublet are expected to avoid the 


tension with BBN induced by the metastable 
charged meson. We left the careful study of BBN 
to future work. 

There are also baryonic composites in 
the SU(N.) confining gauge theory since 
m3(SU(4)?/SU(4)) = Z. Such baryonic compos- 
ites usually saturate the cut off scale 47 fp and 
thus are heavy. They are also the topological 
objects, which in general get more suppressed 
production rate besides of the heavy mass kine- 
matical suppression. Therefore, we only expect 
very tiny production rate at the LHC. In the 
case of My > My2 and Y2= 0 [77], the (~?)N» 
baryonic composites would be a composite DM 
candidate. The thermally averaged annihilation 
rate of the dark baryons could be estimated us- 
ing partial wave unitarity [23, 24, 28], and the 
dark baryon mass could be bounded from above 
as mp < 100 TeV [78]. When this bound is 
saturated, we find fp ~ 2.5 TeV with Nn = 
10 as our benchmark model marked by star in 
Fig. 3. The estimation is made by employing 
the large N scaling for the composite baryon 
mass of mg ~ NnAn [26, 27]. In the early uni- 
verse, once it is thermally produced, the correct 
abundance of baryonic dark matter could be ob- 
tained with a relatively strong coupling. One 
may ask whether this baryonic dark matter is 
metastable through instanton effects [29]. How- 
ever, the enormous suppression factor which is 
proportional to exp(—87?/g?) will make its life 
time much longer than the age of our universe. 
And the consequence of topological dark matter 
on cosmic ray signals and the decay of the DM 
through higher dimension operators is highly re- 
lated with the choice of N, [22, 25], which is be- 
yond the scope of this work. For more generic 
case without considering a hypothetical compos- 
ite DM candidate with Y> = 0 as in the minimal 
model, the 750 GeV diphoton signals can be ac- 
commodated by varying fr from several hundred 
GeV to O(1) TeV with 3 < Nn < 10. 


IV. CONCLUSION 


We have studied the possibility that a sin- 
glet scalar (C'P-even or CP-odd) to account 
for the recent diphoton excess observed by AT- 
LAS and CMS which has attracted a lot of in- 
terests [37, 40-74]. We focused on the gluon- 
gluon initiated process and studied it in an ef- 


fective field theory approach, where the corre- 
sponding higher dimension operators can be gen- 
erated through the heavy vector-like fermions. 
We then consider a natural example that the 
singlet scalar is a pseudo-Goldstone boson from 
chiral symmetry breaking of a new strong sec- 
tor and the interactions with SM gauge bosons 
(the gluon, photon and Z) are purely topologi- 
cal and arising from anomaly. We consider the 
minimal flavor symmetry group SU(4) x SU(4) 
with 15 pNGBs and find the lightest pNGB is 
colour-singlet, which is a good candidate for the 
diphoton excess. Our model also predicts the 
colour-octet, colour-triplet scalars and compos- 
ite baryons. For the colour-triplet, there is no 
topological interaction arising from anomaly due 
to the SU(3), symmetry, therefore it could be 
meta-stable due to higher dimensional operators 
and generates R-hadron like signal at the LHC. 
Other anomaly decays of colour-octet scalars are 
also discussed. The lightest neutral baryon could 
be a viable composite dark matter. If the dipho- 
ton excess is confirmed in the near future, the 
discovery of these resonances and the precise pro- 
duction and decay will provide us a strong test 
on our scenario. 
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Note added. Before this paper is submitted, 
Ref. [37] appeared which also discussed a singlet 
PNGB production and decay through anomaly. 
Nevertheless, our model considers a more general 
charge assignment and focus on the case where 
the second confining fermion 7? is a SM singlet, 
where Nn copies of Y? could be a composite dark 
matter candidate. We also consider the possibil- 
ity that colour triplet scalar is meta-stable, which 
results a R-hadron like signal at the LHC Run 
II (its decay pattern through the higher dimen- 
sional operator is also different). Moreover, we 
calculate the predicted hypercharge of the two 
confining fermions and the confining colour num- 
ber Nn in terms of pNGBs different diboson de- 
cay rate. 


Appendix A: The SU(4) generators and 
WZW effective terms 


We use the generalized Gell-Mann matrices as 
our generators 


1 

1/0 i5 1 0 
Tt == T =e 

(1 o) 2v6 | 0 
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where X* for a = 1,---,8 are the SU(3) Gell- 
Mann matrices with Tr(A*\?) = 28%, and T" 
is the third Cartan generator. T’ with i =9,...14 
are not used in our evaluation and their expres- 
sions are neglected. Collectively, we write down 
the pNGBs as 
IAT = PLT. +T + (PT + DT". 
(A2) 
Our hypercharge generator is defined by 
Y = diag(Y1 , Yı , Yı , Y2). (A3) 


With the conventions listed above, it is 
straightforward to write down the WZW term 
between the pNGB II, and the SM gauge fields 
according to Eq. (7). For the singlet ®,, they 
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where tw = sin w / cos 0w. For octets ®g, they 
read 
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(A5) 


where the symmetric tensor d2°° is given by 


1 
Te Ph = ~ +a% T". (A6) 
The partial decay widths for the singlet ®; can 
be obtained from the WZW term between the 


pNGB II, and gauge fields from Eq. (A4), 
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T[®, > ZZ] = (Y? — Y)’ M3, 


(A7) 


and also the partial widths of the octet ®g from 
Eq. (A5), 


T[®g > gg] = 5 Mès pig > gg] 
8 g= gg 11% + 99), 

Y2 M. 
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